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Abstract. It is reported that a surface model of Polyakov strings undergoes a first-order phase transition 
between smooth and crumpled (or branched polymer) phases. The Hamiltonian of the model contains the 
Gaussian term and a deficit angle term corresponding to the weight of the integration measure dX in the 
partition function. 
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1 Introduction 

A large effort has been devoted to the understanding of the 
phase structure of elastic membrane models of Polyakov- 
Kleinert and Helfrich 3 . The surface models are ex- 
pected to undergo phase transitions between smooth and 
crumpled phases at finite bending rigidity 01|31|SlEIIHlini 
I10iri l,12 13 14 . Numerical studies have also been con- 
centrated on the phase transition llfiUlf 



The bending energy, which is an extrinsic curvature 
term, plays an essential role in smoothing the surface. 
In the ordinary model of membranes, the smooth phase 
emerges owing to the bending energy. Except self-avoiding 
surface models > 25II26II27II28II29I , we have currently no sur- 
face model that has a smooth phase without the extrinsic 
curvature term. 

Baillie- Johnston et. al. |3(JII31II32II33| investigated the 
discretized Polyakov random surface model with intrinsic 
curvature terms in order to find smooth surfaces. One of 
the models was recently investigated in |34| on relatively 
larger lattices, however, there was no smooth phase in the 
model, although branched polymer like phase was found. 
We consider that the question about a smooth phase in the 
models remains unanswered, and the surface model with 
intrinsic curvatures remains to be studied. Thus, moti- 
vated by the attempts of Baillie-Jhonston et. al., we study 
in this article a model with an intrinsic curvature term and 
see that the smooth phase and the crumpled phase are ap- 
peared and connected by a phase transition. Our results 
may lead to better understanding of the role of curva- 
ture energy in surface models undergoing phase transi- 
tions and/or shape transformation, although the surfaces 
are allowed to self-intersect and hence phantom. 

Two kinds of models are investigated; the tethered 
model and the fluid model. The tethered model is defined 



68.60.-p Physical properties of thin films, nonelec- 



on fixed connectivity surfaces, and the fluid model on dy- 
namically triangulated surfaces. The Hamiltonians of the 
models are identical with each other. It will be shown 
that the tethered model has the smooth and the crum- 
pled phases, which are connected by a first-order phase 
transition. Moreover, we will find in the fluid model that 
there appears a new phase, which resembles the branched 
polymer phase |35II36| . between the smooth phase and the 
crumpled one. We understand from the results obtained 
in this article that the smooth phase is stable against such 
as thermal fluctuations in each model defined only with 
intrinsic variables of the surface. 



2 Model 

The partition function of the fluid model of fixed number 
of vertices is defined by 



N 



(1) 



T " i=l 



5(X,r) = 5i-a^3, 



where the center of the surface is fixed to remove the trans- 
lational zero mode. S(X^ T) denotes that the Hamiltonian 
S depends on the position variables X of vertices and the 
degree of freedom for the triangulation T. S\ and 6*3 are 
defined by 

5i=^(X,-X,f, 53 = ^log(5,/27r), (2) 

where i'^ "^i svcav over bonds (if) connect- 

ing the vertices i and j. The bonds (ij) are edges of 
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the triangles. Si in Eq. (jSJ is the vertex angle, which is 
the sum of the angles meeting at the vertex i. Recall- 
ing that 5i — 2ti is the deficit angle, we call the deficit 
angle term. We note that the sum of the deficit angles 
— 27r) is constant on surfaces of fixed genus be- 
cause of the Gauss-Bonnet theorem. The partition func- 
tion of the tethered model is, on the other hand, defined by 
Z = j n»=i dX, exp [-S{X, T)], where T is fixed. Both Si 
and S3, reflect intrinsic properties of the surface. Note also 
that Si/N = 3/2, which comes from the scale invariance 
property of Z in both models. 

We should refer to a relation between the deficit an- 
gle term S^ and the integration measure dXi of the parti- 
tion function in Eq. From a conformal field theoretical 
viewpoint P7ll38ll39llin| , the integration measure Y\^ dXi 
can be replaced by qfdXi, where qi is the co-ordination 
number of the vertex i. This a is believed to be 2a — i. 
On the other hand, qf is considered as a volume weight 
of the vertex i in the measure dXi. Hence, it is possible 
to extend 2a to non-integer numbers by assuming that 
the weight can be chosen arbitrarily. Moreover, it is also 
possible to extend qi to continuous numbers because of 
the same reason. Hence, the weight qf can be replaced 
by Oi^fi which can also be written as exp(a logJi). 
Including a constant weight (27r)", we have 5*3 in Eq. 

The model in this paper is very similar to the one 
investigated in [30113 1II32II33II34| . where the fluid and the 
tethered models were interpolated. A model in 30 ,31,33 
15505^ includes the term a^ - \qi — 6| or a^^{qi — 6)^. 
Using the expression log((5i/27r) = log[l-|-((5i — 27r)/27r] and 
replacing (5^ by qi and 27r by 6 in S^, of Eq. 121) , we have the 
term proportional to ~ + (higher order). Thus, 

we confirmed that the lowest order term in Sj, is identical 
with the expression Q;^j(gi — 6)^ described above upto 
multiplicative constant. We note that the similarity be- 
tween 5*3 of Eq. (O and a deficit angle term denoted by 
•S'tight in Pni can also be seen. 

We will see below that the fluid surfaces become crum- 
pled, branched polymer like, and smooth, when a increases 
from a = 0. A higher-order transition separates the branched 
polymer phase from the crumpled one, and a first-order 
transition separates the branched polymer phase from the 
smooth one. On the other hand, we find in ^31^ that the 
fiuid surfaces become crumpled, branched polymer, and 
crumpled, with increasing a. One of the transitions is of 
higher (or second) order, and the other is of first order, in 
the model of 

It should also be noted on differences between 5*3 (5) of 
Eq. ISJ and the corresponding term S'3(g) = log(gi/6), 
which can be expanded to X]j(9j ~ 6)^ -|- (higher order) as 
described above. Although the difference between these 
two S-i comes only from the one between q and (5, the 
roles of q and 6 in the model make large differences in 
the results. A major difference between the fluid model 
in this paper and the one in [35 is whether the model 
has a smooth phase or not. q is unchanged in the MC 
update of X, whereas 5 changes in that process. As a con- 
sequence, while the co-ordination dependent term 53 (g) 
is non-trivial only in the fluid model, S^iP) in Eq. is 



non-trivial both in the tethered and in the fluid models. In 
fact, the model with X]j('?j ~6)^ has no smooth phase |34j . 
whereas the model with 6*3 ((5) in Eq. jSJ has the smooth 
phase both in the tethered and in the fluid surfaces as we 
will see below. 



3 Monte Carlo Simulations 

The spherical surfaces, the starting configurations of MC 
for the fluid model, arc discretized by the Voronoi trian- 
gulation in R'^. The radii of the initial spheres are 
chosen so that S'i/iV~3/2. The surfaces are constructed 
relatively uniform in co-ordination number; the maximum 
co-ordination number is about 8. 

On the contrary, the uniform lattices, on which the 
tethered model is deflned, were obtained by MC with dy- 
namical triangulation for a model defined by 5 = S'i-l-55'2— 
aX)i log(?i/6), where Si = X]j(l ~ cos6'i) is the bending 
energy. Using the parameters a = 280, 6 = 2 (or a = 350, 
& = 2), we obtained the uniform lattices. By monitoring the 
number of vertices of = 5, we stopped the MC run 
when first reduced to 12, and thus obtained the uni- 
form lattices. The uniform lattices, which will be utilized 
in simulations for the tethered model, are thus character- 
ized by iV5 = 12 and A'6 = 3iV — 18. In each lattice, only 12 
vertices are of qi = 5, and all other vertices are of qi = 6. 
Except the lattices of A^s = 12, there exists no lattice whose 
qi are of = 5 or = 6. This is easily understood from the 
fact that A^B = 3A^— 6 on triangulated surfaces of spherical 
topology, where A'b is the total number of bonds. 

The canonical Metropolis MC technique is used to up- 
date the variables X and T. In the MC simulations, the 
new position of X is chosen so that X' =X + 5X, where 
5X takes value randomly in a small sphere. The radius 5r 
of the small sphere is defined by using a constant number 
e as an input parameter so that 

5r = e{l), (3) 

where {I) is the mean bond length computed at every 250 
MCS (Monte Carlo sweeps), and {I) is constant in the 
equilibrium configurations because Si/N = 3/2. The new 
position X' is accepted with the probability 
Min[l,exp(-Z\S')], where ZiS* = 5(new) -5(old). The pa- 
rameter e in Eq. (PJ is chosen to maintain about 50% ~ 
60% acceptance for X. A lower bound 10~^ x Aq is as- 
sumed for the area of triangles, where Aq is the mean 
area of triangles evaluated at every 250 MCS. X' is gen- 
erated so that the resulting areas of triangles are larger 
than IO^^xAq. It should be emphasized that the results 
obtained by using Sr in Eq. JSJ arc identical with those 
by using a constant Sr. 

The triangulation T is updated by the bond fiip. Then, 
the acceptance rate rq- for the bond flip is uncontrollable. 
We have, in our MC simulations, 15% < rq- < 35%, and 
flnd that rq- discontinuously changes at a discontinuous 
transition which will be shown below. 

The bond is flipped so that the resulting areas of trian- 
gles are larger than 10~^x as in the process for X. No 
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restriction is imposed on the bond length in the update 
of the variables X and T. The maximum co-ordination 
number q™^^ of the fluid surface of size N — 1500 was 
about grmax ^ 20 both in the smooth phase and in the 
branched polymer phase, and q'^'^'^ ~ 30 in the crumpled 
phase. The histogram h{q) of the co-ordination number q 
will be shown below. iV-trials for X and iV-trials for T 
make one MCS in the fluid model. 

The canonical ensemble average {Q) is computed via 
the time series average of a series {Qi}. A series {Qi} is 
sampled at every 500 MCS after thermalization sweeps, 
which are large enough so that the mean square size X^, 
defined by 



X' 



X) 



X 



(4) 



becomes statistically independent from the starting con- 
figurations. In fact, 0.5 x 10'' MCS are discarded for the 
thermalization of surfaces of iV < 600, and 1 x 10^ MCS 
for N > 1000 surfaces. Total number of MCS is about 
(8 ~ 10) X 10^ at the vicinity of the transition point of 
N > 1000 fluid model and of TV > 1500 tethered model, 
and it is relatively smaller in other cases. 



The rate of acceptance rx for X and rx for bond flip 
are plotted in Fig.^b). We see in Fig.Q^b) that rx is al- 
most constant and rx is discontinuous. In order to main- 
tain 0.5 < rx < 0.6, we adjust e in Eq. to a suitable 
value, because rx rapidly changes against a at the tran- 
sition point. This is the reason why rx is not completely 
constant. The discontinuous behavior of rx also indicates 
that the model undergoes a discontinuous transition. 
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Fig. 2. The normalized histograms h{q) of the co-ordination 
number q of the fluid surface of A'^ = 1500 at (a) q = 10 (crum- 
pled phase), (b) a = 500 (branched polymer phase), and (c) 
a =1600 (smooth phase). 
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Fig. 1. (a) AiX^) at a = 1400 and at a = 1600 on = 1500 
fluid surface, and (b) the acceptance rates rx for X and tt for 
T of the fluid surface model. 



The model membrane is expected to swell out at large 
Q. and crumple at small a. Therefore, the shape of mem- 
brane, which is an extrinsic property of surfaces, is charac- 
terized by X"^ . In order to see the correlation in the series 
of X^, the autocorrelation coefficient 



E,X2(r,)X2(r, 



A(X2 

-1 = Ti + rt X 500, n = 1, 2, ■ 



(5) 



is plotted in Fig. QJa), where X'^iji) is sampled at every 
n X 500 (n = 1, 2, ■ • ■) MCS. We easily see in Fig.IJa) that 
A(X^) at a = 1400 and the one at a = 1600 are completely 
different from each other. This result implies that the con- 
vergence speed of our MC on = 1500 fluid surfaces at 
a = 1400 is about 10 times slower than that at a= 1600, 
and indicates the existence of some discontinuous phase 
transition. 



In order to see the distribution of the co-ordination 
number the normalized histograms hijq) are shown in 
Figs.|ll;a),[2l;b), andlUc). The histograms /i(g) in the fig- 
ures were obtained in the final 5 x 10^ MCS on the fluid 
surface oi N = 1500 at (a) a = 10, (b) a = 500, and (c) 
a = 1600. These three values of a correspond to three 
distinct phases; the crumpled, the branched polymer, and 
the smooth phases, which will be clarified later. We note 
that the co-ordination number q does not always reach the 
above mentioned qf^^'^ the maximum co-ordination num- 
ber during the final 5 x 10^ MCS in each case. Thus, we 
find from Figs.[21[a),|21[b), andEfc) that the co-ordination 
number q is dominated by 4 < g < 7 at each a in the fluid 
model. No distinct difference can be seen in h{q) obtained 
in the three phases, although it can bee seen that h{q = 3) 
in the crumpled phase at a = 10 is larger than those in the 
other two phases. It is remarkable that h{q) in the smooth 
phase shown in Fig. EIc) is almost identical with h(q) in 
the branched polymer phase shown in Fig. ^[h) , because 
we expect that the curvature energy 5*3 is closely related 
with the co-ordination number q and plays an essential 
role in the phase transition. The variation of q is also very 
small even at the discontinuous transition point between 
the branched polymer and the smooth phases. 

Nevertheless, we can find in S3/N a discontinuous change, 
which reflects a discontinuous transition between the branched 
polymer and the smooth phases, in both models. However, 
the gap in S3/N is very small and comparable with the 
standard deviation of S3/N. Hence, we can not conclude 
only from the gap that the model undergoes a discontin- 
uous transition. 

Si/N vs. b are plotted in Figs.Ofa) andOl^b). We see 
in Fig.|2Ia) that the expected relation Si/N = 3/2 is sat- 
isfied in the tethered model, and also in Fig. Olb) that 
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Fig. 3. (a) Si/N obtained from = 1500 and iV = 2500 
tethered surfaces, and (b) Si/N obtained from A'^ = 1000 and 
AT = 1500 fluid surfaces. 



Si/N — 3/2 in the fluid model. The discontinuous change 
seen in Si/N of TV = 1500 in Fig. Olb) implies the exis- 
tence of a discontinuous transition. Since the gap is almost 
equal to or less than 0.4% of Si/N, we consider that the 
relation Si/N — 3/2 is not influenced by the discontinu- 
ous transition. Thus, we understand that the surfaces in 
both sides of the transition point are in the equilibrium 
configurations prescribed by Z of Eq. . 
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Fig. 4. vs Q of (a) tethered surface model and of (b) fluid 
surface model. Dashed lines in (b) denote the phase boundaries 
on 7V= 1500 fluid surfaces. 



The mean square size X'^ is plotted in Figs. EJa) and 
^b). In the tethered model, X'^ is obtained on surfaces of 
size up to A^ = 2500. We see a gap in each in the teth- 
ered model. These discontinuous changes of imply that 
the phase transition is of first order. On the other hand, 
the dependence of X^ on a in the fiuid model appears to 
be continuous on surfaces of iV=1000 and iV=1500. The 
reason why X^ takes intermediate values at a just below 
the transition point is that the surface belongs to a new 
phase, which emerges only in the fluid model. This new 
phase resembles the branched polymer phase pnii:i5li;^ 6' in 
shape of the surfaces, which will be shown as a snapshot 
below. Thus, we call the phase as the branched polymer 
phase, although we can hardly conflrm that the new phase 
is exactly identical with a branched polymer phase from 
the present numerical results. 



It should also be noted that the fluid surfaces become 
crumpled when a — > 0, as seen in Fig. E^b). The dashed 
lines drawn vertically in Fig.^b) denote the phase bound- 
aries of the fluid surfaces of size TV = 1500. The leftmost 
data X^ in Fig. Efb) were obtained at a = 10 on surfaces 
of iV = 1000 and N = 1500, and their dependence on N 
clearly different from those obtained at larger a. Thus, we 
considered that the phase boundary between the crum- 
pled and the branched polymer phases is located close at 
a = lG. 
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vs TV on (a) tethered surfaces and on (b) fluid 
surfaces. denoted by smooth in (a) and (b) were obtained in 
the smooth phase close at the transition point in both models. 
X^ denoted by crumpled in (a) were obtained just below the 
transition point. X^ of surfaces < 600 denoted by crumpled 
in (b) were obtained just below the transition point, and those 
of Af > 1000 were obtained at q = 10. 



Figures lUa) and[Sfb) are log-log plots of X^ against 
N. X^ denoted by smooth in the figures were obtained 
in the smooth phase close at the transition point in both 
models. X^ denoted by crumpled in the tethered model 
were obtained just below the transition point. X^ denoted 
by crumpled in the fluid model of A'^ < 600 were obtained 
just below the transition point, and those of A^ > 1000 were 
obtained at a = 10. The straight lines give the Hausdorff 
dimension H, which is defined by 



X^ NTT, 



(6) 



The surface swells and becomes almost smooth in the 
smooth phase both in the tethered model and in the fluid 
model; we already confirmed this from the data in Figs. 
0Ja) and^fb). X^ has a value of the order of the radius 
squares when the membrane becomes a swollen sphere. 
Hence, we expect that iJ 2 in the smooth phase. We 
have in fact 



H = 2.11± 0.03 (tethered, smooth), 

H = 10.47 ± 1.23 (tethered, crumpled), 

H = 1.99 ± 0.07 (fluid, smooth), 

i/ = 4.15 ±0.47 (fluid, crumpled). 



(7) 



The results in Eq. iQ obtained in the tethered model 
indicate that the phase transition between the crumpled 
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and the smooth phases is characterized also by a gap of H. 
The reason why H is finite in the crumpled phase comes 
from the fact that X"^ are obtained close at the transition 
point, which is far distant from a = 0, where H is expected 
to be infinite. 

In the fluid model, on the contrary, the dependence 
of on a in the branched polymer phase can not be 
extracted from the present numerical data. In order to 
see how H depends on a in the branched polymer phase, 
large-scale simulations are necessary. Nevertheless, it is 
possible that H continuously changes against a in the fluid 
model, if the new phase is the branched polymer one. In 
fact, it is expected that ~ 2 in the branched polymer 
phase [24,35 ,36 . In addition, we see in Fig.Q^b) that 
smoothly changes between the branched polymer phase 
and the crumpled phase. 
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Fig. 6. S2/NB vs b obtained on (a) tethered surfaces and on 
(b) fluid surfaces. S2 is defined by ^'2 = y^^(l — cos Oj). 



Figures El^a) andE^b) are the bending energy S2/NB 
vs a respectively obtained on the tethered surfaces and 
on the fluid surfaces, where S2 is defined by S2 — J2ii^~ 
cos 9i), Nb the total number of bonds. 5*2, which is an 
extrinsic variable and is not included in the Hamiltonian, 
represents how smooth the surface is. The gap seen in 
S2/NB also indicates that the phase transition is of first 
order. Thus, we see a first-order transition between the 
smooth and the crumpled phase in the tethered model, 
and also see a first-order transition between the smooth 
phase and the branched polymer phase in the fluid model. 

We note that first-order nature of transitions is ordi- 
narily confirmed by a discontinuous change (or a gap) of 
energy term included in Hamiltonian. A gap in some phys- 
ical quantity that is not included in Hamiltonian can also 
confirm loosely a first-order nature of transitions. Thus, 
we consider that the gap in in Figs. Ela) and^fb), 
and the gap in 5*2 in Figs. EJa) and^Jb) indicate the ex- 
istence of first order transitions. 

Figures [7|[a) , [7fb) and [7fc) are snapshots of surfaces 
obtained respectively at the smooth phase, at the branched 
polymer phase, and at the crumpled phase in the fluid 
surface model. The tethered surface model gives the same 
appearance of surfaces as the fluid surface model both in 




1(c) Of =1600 




Fig. 7. Snapshots of fluid surface obtained in (a) the crumpled 
phase at a = 10, (b) the branched polymer phase at a = 1200, 
and (c) the smooth phase at a= 1600. Small spheres represent 
the vertices. Af = 1500. 



smooth and in crumpled phases, although the tethered 
model has no branched polymer phase. 



4 Summary and Conclusions 

We have studied extrinsic properties of surface models of 
membranes embedded in R^. The Hamiltonian contains 
the Gaussian energy term 6*1 and the deficit angle term 
S3; S = Si — aSs- This 5*3 was introduced by a straight- 
forward extension of the co-ordination dependent term, 
which comes from the weight of the integration measure 
dX. Two kinds of models were investigated: one is a teth- 
ered model and the other is a fluid one. We found in the 
tethered model that there are two distinct phases; smooth 
and crumpled, and that these phases are connected by 
a discontinuous transition. It was also found in the fluid 
model that there is a new phase (branched polymer phase) 
between the smooth and the crumpled phases. A discon- 
tinuous transition separates the branched polymer phase 
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from the smooth one in the fluid model. Although we 
called the new phase in the fluid model as the branched 
polymer phase, large scale simulations are necessary to 
confirm that the new phase is exactly identical with the 
branched polymer one. 

This work is supported in part by a Grant-in- Aid for 
Scientific Research, No. 15560160. 
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